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ABSTRACT 
The Lorentz factor (LF) of gamma-ray burst (GRB) ejecta may be constrained by observations of high-energy 
(HE) spectral attenuation. The recent Fermi-LAT observations of prompt GeV emission from several bright 
GRBs have leaded to conclusions of unexpectedly large LFs, T > 10 3 . Here we revisit this problem with two 
main concerns. (1) With one-zone assumption where all photons are assumed to be generated in the same region 
(radius) and time, we self-consistently calculate the 77 optical depth by adopting a target photon spectrum with 
HE cutoff. We find that this might be important when the GRB LF is below a few hundreds. (2) Recent Fermi- 
LAT observations suggest that the bulk MeV-range and HE (> 100 MeV) emission may arise from different 
regions. We then consider a two-zone case where HE emission is generated in much larger radii than that of the 
MeV-range emission. We find that the HE emission may be mainly attenuated by MeV-range emission and that 
the attenuated HE spectrum does not show an exponential spectral cutoff but a slight steepening. This suggests 
that there may be no abrupt cutoff due to 77 attenuation if relaxing the one-zone assumption. By studying the 
spectra of three bright Fermi-LAT GRBs 0809 16C, 090510 and 090902B, we show that a bulk LF of T ~ 600 
can be consistent with observations in the two-zone case. Even lower LFs can be obtained in the multi-zone 
case. 
Subject headings: Gamma ray burst: general 



1. INTRODUCTION 

Relativistic expansion is a key property of gamma-ray 
bursts (GRBs), and has been confirmed by measurements of 
radio afterglow sizes, for examples, the indirect estimation by 
radio scintillation in GRB 970508 (Waxman et al. 1998) and 
direct imaging of nearby GRB 030329 (Taylor et al. 2004). 
These observations revealed mildly relativistic GRB ejecta, 
r ^a few, in the radio afterglow phase. However, it is well be- 
lieved that GRB ejecta are ultra-relativistic in the beginning- 
this is required to solve the so-called "compactness problem" 
(e.g., Piran 1999). The compact GRB source, suggested by 
the rapid variabilities in MeV light curves, and the huge lumi- 
nosity suggest hot, optically thick GRB sources, which is in 
confliction with the nonthermal and hard GRB spectra. Rela- 
tivistic expansion of the emission region is introduced to solve 
this problem. In order for the ^100 MeV photons, as detected 
by EGRET in several GRBs, to escape from the emission re- 
gion, avoiding 77 attenuation, the bulk Lorentz factor (LF) of 
the emission region is required to be extremely large, T > 10 2 
(e.g., Lithwick & Sari 2001; Krolik & Pier 1991; Fenimore et 
al. 1993; Woods & Loeb 1995; Baring & Harding 1997). Re- 
cently, the powerful Fermi satellite reveals in much more de- 
tail the high-energy (HE) emission from GRBs. Several bright 
GRBs are reported to show time-integrated spectra extending 
up to GeV or even tens GeV, without any signs of spectral 
cutoff. Assuming the 77 optical depth for these HE photons 
are below unity, these observations have leaded to even larger 
bulk LFs, T > 10 3 (Abdo et al. 2009a,b,c). This is putting the 
theoretical problem of relativistic jet formation to extremes. 
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In the previous constraints two assumptions are usually 
taken. First, all photons, from low to high energy, are pro- 
duced in the same region and the same time. This "one-zone" 
assumption is not solid, as Fermi observations actually re- 
vealed that: the onset of HE emission is delayed relative to 
MeV emission (e.g., Abdo et al. 2009a,b,c); the HE emission 
lasts longer than MeV emission (e.g., Abdo et al. 2009a,b,c); 
the bulk emission shifts toward later time as the photon en- 
ergy increases (Abdo et al. 2009a) and the shift is longer than 
the variability times in MeV light curves, as pointed out by 
(Li 2010); some GRBs obviously show distinct HE compo- 
nents with different temporal behaviors (Abdo et al. 2009b, c). 
All these features may imply that different energy photons are 
produced in different regions. 

In particular, the bulk > 100 MeV emission in GRB 
0809 16C shows ~1 s shifting relative to MeV emission, 
which is much longer than the MeV variability time, < 100 ms 
as revealed by INTEGRAL (Greiner et al. 2009), strongly im- 
plying that > 100 MeV emission is produced in a region of 
much larger radii than MeV emission's (Li 2010). As pointed 
out by Li & Waxman (2008), within the framework of inter- 
nal shock model, the internal collisions at small radii, which 
would produce the prompt MeV emission, are expected to 
lead to "residual" collisions at much larger radii, which would 
produce low-frequency emission. The electrons accelerated 
by residual collisions at larger radii inverse-Compton scat- 
tering the MeV photons and/or double scattering the low- 
frequency photons could produce HE emission (Li 2010; 
Zhao et al. 2010). In this case, MeV and HE photons are 
produced in different regions. In the comoving frame of HE 
emission region, the MeV photons would be collimated other 
than isotropic, thus the 77 absorption is angular dependent. 

Second, the target photon spectrum is assumed to be ex- 
tending to infinity. As pointed out by Li (2010), the calcu- 
lation of 77 optical depth taking such a target photon filed 
is obviously not self-consistent, because the HE spectral end 
should be cut off due to absorption considered in the calcula- 
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tion. 

In this paper, we revisit the problem of GRB LF constraint 
by modifying the above mentioned two assumptions. We con- 
sider in §2 a one-zone case where the 77 optical depth is cal- 
culated self-consistently by assuming a truncated target spec- 
trum, then we consider in §3 a simple two-zone case with 
anisotropic effect on 77 optical depth taken into account. In 
§4 we studied the spectra of the three bright Fermi-LAT GRBs 
and constrain their LFs. §5 is discussion and conclusions. 
In the following we assume the concordance universe model 
with (n m ,n A ) = (0.27,0.73) andtf = 71kms _1 Mpc- 1 . 

2. ONE-ZONE CASE 

Consider a GRB ejecta with bulk LF T and radius R. As- 
sume the photons in the comoving frame of the ejecta is 
isotropic, with photon number density per photon energy 
dn' /de' . Hereafter, unless specified otherwise, quantities with 
prime denote the comoving frame, and non-primed ones de- 
note the frame of observer on the Earth. 

In the (comoving-frame) dynamical time R/Tc, a photon 
travels a path of R/T, For a photon of energy e' = e(l +z)/T 
(with z the GRB redshift), the optical depth due to 77 colli- 
sions during a dynamical time is given by (Gould & Schreder 
1967) 
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where jj,' = cos 6' and 9' is the angle between the colliding 
photon pair. The cross section is given by 
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where j3 e = ^ I - (m e c 2 / E) 2 and E = y/e'e'(l—ji')/2 are the 
velocity and energy, respectively, of the generated electron in 
the center of momentum frame of the collision. The radius R 
of the emission region can be related to the angular spreading 
time <5? ang , due to geometry effect, by R = 2r 2 c<5f ang /(l+z). 
As the angular spreading time is related to the observed vari- 
ability time St by 5t ang = St, we have 
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For a GRB with the observed photon number per unit time per 
unit photon energy per unit detector area, denoted by N(e), the 
photon number density per unit photon energy in the comov- 
ing frame can be given by 
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where dh is the GRB luminosity distance, and e = IV/(1 +z). 
It is important to note a difference from the previous works. 
In eq. (1) we did not take the upper limit of the integration to 
be infinity but a certain photon energy e' miix , because the HE 
tail is expected to be cut off due to 77 absorption. The cutoff 
energy is just where T(e' miix ) = 1 happens. To self-consistently 
solve out the cutoff energy e max = Te' m , dx /(l +z) for given T, 
we need to take the upper limit of the integration to be e^ ax , 
and solve r(ej nax ) = 1 using eqs. (1-4) and observed GRB 
spectrum N(e). 



It is well known that the GRB spectrum can be fit by the 
Band function (Band et al. 1993) 
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where e c = e p (a — /3)/(2 + a), and A, a, j3 and e p are the nor- 
malized coefficient, low-energy slope, HE slope and the vF v 
peak energy, respectively. In some Fermi-LAT GRBs an extra 
spectral component beyond the Band-function is claimed to 
exist, especially in HE end (Abdo et al. 2009b, c). This extra 
component can be described as a power law, 
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with Apl the normalization at 1 GeV and /3pl the spectral in- 
dex. 

It is helpful to solve out the r-e max relation with some 
approximations first. Typically the HE, > 100 MeV, photons 
mainly interact with photons above the peak energy. Let us 
approximate the target photon distribution as a single power 
law N(e) = Noe~ s in the following analytical derivation. 

In eq. ( 1 ), usually the upper limit of the first integral is taken 
to be 00. This is valid for e max 3> r 2 m 2 c 4 /[e max (l +z) 2 ] and 
the spectrum slope s > 1 . In this case, using ^-approximation 
for the cross section at target photon energy above the thresh- 
old, a fa (3/16)0-7-, r(e max ) = 1 can be solved to give T as 
function of £ max , 
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However, when e max > r 2 m 2 c 4 /[£ max (l +z) 2 ], i.e., the energy 
of annihilated photons is compared with that of target pho- 
tons, the upper limit cannot be taken as 00 any more. In the 
case, r is given by (Li 2010) 



rd+z). 
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Next we carry numerical calculation to solve out T(e m3x ) = 
1 . For the observations, we take the three bright Fermi-LAT 
GRBs 080916C, 090510 and 090902B, and consider the same 
time intervals in the GRBs where the LFs have been con- 
strained by Abdo et al. (2009a,b,c), as well as section a in 
GRB 080916C. The properties of spectra and flux for these 
GRBs are shown in Table 1. The calculated results are given 
in Fig 1 , where we compare the results of self-consistent cal- 
culation and previous method using a target photon spectrum 
without HE cutoff. We see that the results deviate each other 
for e max < 100 MeV or T <a few hundreds. In the case of sec- 
tion a in GRB 0809 16C, where the maximum observed pho- 
ton energy is lower (see Fig 1), the Lorentz factor limit with 
the self-consistent calculation is much smaller than that with 
the previous method. Thus to be self consistent, the upper 
bound of the integration in eq (1) should be carefully taken as 
the maximum photon energy in this case. We also note that 
the LF constraints using upper limit of infinity are still valid 
for those time segments that have been used by Abdo et al. 
(2009a,b,c). 

3. TWO-ZONE CASE 

As discussed in the introduction, the Fermi-LAT observa- 
tions hint that there may be different emission regions of dif- 
ferent radii in GRB prompt emission. As the HE delay of 
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TABLE 1 

The parameters of three bright LAT-GRBs 



GRB name 


Time interval 


e i> 


a 





A 


A>L 


Ap L 


z 


8t 


^highest 




(s) 


(keV) 






(cnrV'keV- 1 ) 




(cnrV^eV- 1 ) 




(ms) 


(GeV) 


GRB 080916C-a 


0.004-3.58 


440 


-0.58 


-2.63 


0.055 


_ 


_ 


4.35 


100" 


0.02* 


GRB 080916C-b 


3.58-7.68 


1170 


-1.02 


-2.21 


0.035 


- 


- 


4.35 


100° 


3 


GRB 090510 


0.8-0.9 


1894 


-0.86 


-3.09 


0.028 


-1.54 


6.439 x 10-' 


0.903 


12 


30.5 


GRB 090902B 


9.6-13 


821 


-0.26 


-5.0 


0.082 


-1.98 


4.3xl0- 10r 


1.822 


53 


11.2 



NOTE. — References: a: Greiner et al. (2009); b: the spectral flux of this section at > 20 MeV is only upper limit, as shown in the supporting material of 
Abdo et al. (2009a); c: private communication with Francesco de Palma; and the other parameters are taken from Abdo et al. (2009a,b,c). 

comoving frame of HE emission region, the solid angle is then 
Afi' = 27r(l-cosa'X with cose/ = (cosa-/?r)/(l _ /3rcosa) 

and /3 r = VT 2 -l/T. 

The optical depth is not only energy-dependent but also 
angle-dependent. Consider a HE photon of e' travelling with 
an angle 9' (p! = cos 9') with respect to the central axis of the 
target photon beam, then the optical depth corresponding to 
the distance it travels in a dynamical time is given by 




1-GRB080916C(b section) 
2-GRB 080916C (a section) 
3-GRB 090902B 
4-GRB 090510 

- - non-cutoff target spectrum 

— self-consistent 



10 B 
(keV) 



FIG. 1 . — The relation between the observed maximum photon energy and 
the lower limit to the bulk LF in the one-zone case for the three bright GRBs. 
The adopted parameters of the GRBs are shown in Table 1 . As marked in the 
plot, the dash lines correspond to results using target photon without spec- 
tral cutoff, while the solid lines correspond to our self-consistent calculations 
using truncated target photon spectra. The stars denote the observed highest 
energy of photons in the relevant time intervals. 

onset and the shifting of the bulk HE emission are in seconds 
scale, whereas the MeV-range variability times, reflecting the 
dynamical time of the MeV emission region, are in tens of 
ms scale, the MeV emission regions have much smaller, by 
orders of magnitude, size (radius) than that of HE emission 
regions. 

Consider that the ejecta expand to radius R where HE emis- 
sion is being produced. The photons that are emitted in much 
smaller radii and just arrive at radius R should be produced by 
those ejecta released from the central engine with a time delay 
of td = R(l +z)/2T 2 c. If td > St then td is also the observed 
delayed time scale of the HE emission. Thus, once we ob- 
served a time delay tdC^> St) for HE emission relative to MeV 
emission, the HE emission size is implied to be 



R = 2T 2 c 



td 
\+z 



(9) 



As the MeV emission comes from inner regions with 
smaller radii, /?MeV = 2r^ eV c(5f /(l +z) <C R, the MeV photons 
in the comoving frame of trie HE emission region are beamed. 
Here we also denote the LF of MeV emission region as TMeV 
since it may be different from the one of the HE emission re- 
gion, T, in the framework of internal shock model, and the 
difference could be small, TMeV ~ T. 

Consider the geometry plotted in Fig 2. Due to the relativis- 
tic beaming effect, the MeV emission beam that illuminating 
a HE photon produced at R can be approximated as a "MeV 
photon cone" with half open angle of a = Rmsv/R^m^v- Out- 
side of the cone the MeV photon flux can be neglected. In the 
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Here ft' = cos 0' with 0' the angle between HE photon and 
the colliding target photon, d 2 n' /de'dfl' is the energy dis- 
tribution of target photons per unit solid angle, with dCl' = 
sin O'dO'dqy', and <f>' = -k if 6' < a' and ft > cos(a' - $'), oth- 
erwise 



: arccos 



cos a' — cos 0' cos 9' 
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The comoving frame target photon density dn' /de' is given 
by eq (4) but with e = 2IV/(1 +z)- the factor 2 appears for 
the highly beamed case of a' <C 1 . 

In the extreme case when RMev/R —> or TMeV/T —¥ +oo 
the target photons are totally beamed (hereafter other cases 
are called partly beamed) in the comoving frame of the HE 
emission region, then Ail' —J- and the optical depth reduces 
to 

r(eV) = ^/ de'^(l-n'ME). (12) 

1 J2m 2 c 4 /s'(l-fJ,') " e 

Consider an area element in the sphere emitting photons 
which lies at an angle 9 with respect to the line of sight, then in 
its comoving frame a photon travelling along line of sight has 
an angle with respect to the central axis of the target photon 
beam of 

" = T^r, (13) 

Due to Doppler effect, the photon energy in the comoving 
frame is related to the observed photon energy as 

e' = T(l-p r ^e(l+z). (14) 

Denote d 3 P' /dVl'de'dS as the (comoving-frame) emitting 
power per unit solid angle per unit photon energy by material 
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MeV emission radius 



HE emission radius 



In the comoving frame 
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MeV photon axis 




6'>a' 



FIG. 2. — The schematic diagram of the geometry when the a HE photon 
produced at large radius R is being illuminated by the MeV photon beam from 
much smaller radii RueV ■ The upper panel is for the observer frame, while 
the bottom two panels for the comoving frame of GeV emission region, with 
8' < a' and 6' > a' respectively. 

in per unit area of the sphere surface. This emission should 
be modified by 77 attenuation factor e" T(e •>* \ The observed 
"time-averaged" flux is, then, integration over the sphere, 



F e oc / dS 
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Assume isotropic emission power in the comoving frame, 
constant emissivity along the sphere, and power law depen- 
dent on photon energy, then 
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Using dS = 2irR d[i and eqs. (9),(13) and (14), we have 



F e (xtJT- h+ h- h+1 f(e;T), (17) 



where 



/(e;r)= di41-PriiT- l e 



h-\ -T(e',//) 
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Note that /(e; T) is the suppression factor of the primary spec- 
trum. 

It is useful to analyze this factor analytically for the totally 
beamed case. As shown in Appendix, for the single power- 
law target photon distribution, N(e) = Noe~ s and using approx- 
imation a{E) w <jqE~ 2 , the / factor can be approximated as 
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with the break energy at 
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We carry numerical calculation of / factor, and show the re- 
sult in Fig 3. The analytical result is a good approximation. 

Thus in the totally beamed case the spectrum is not affected 
until e > £b,, where the spectrum steepens by a factor of - - 
1. Thus unlike in the case of isotropic target photons, the 
spectrum is not cut off exponentially but show a steepening 
power law. This can be easily understood- in the beam target 
photon case, the HE photons always can escape if they travel 
with a small enough angle with respect to the target beam. 

The break energy is LF-dependent, thus the detection of the 
break in the spectrum can be used to measure the LF of GRBs. 
For GRB 0905 10, the break energy is £ br « 1 GeV for T = 600 
and td = 0.1s. 

In Fig 3 we also show the numerical results for partly 
beamed cases. As can be seen, for the energy range of inter- 
ests, say, < 1 TeV, the partly beamed cases with Tu e \ = T ap- 
proaches the totally beamed case when Rmzv/R < 0.1, which 
is just the case we are considering because St <C fj. We 
also illustrate the small effect of LF variation by showing the 
cases of rMev/r = 0.5-2 with RMev/R = 0.01. Indeed when 
Tiviev > r as expected, the MeV photons are more strongly 
beamed then the situation is more approaching the totally 
beamed case. Even when TMev/r < 1 there is only very little 
effect at very high energy (see Fig 3) since RMev/R "C 1 and 
the MeV photons are still highly beamed. Thus we will ignore 
the effect of variation of LFs and only consider TmcV = T in 
the following calculations. 

It should be noted here that in the two-zone case, besides 
the HE absorption due to the inner-coming beamed MeV pho- 
tons, the absorption due to interactions with photons locally 
originated from the HE emission region can also contribute 
to the total optical depth. This adds an extra attenuation fac- 
tor e" T »if( £ ) in the resulted spectrum, where the optical depth 
T se if(e) is given by eq (1) with R being the HE emission region 
radius eq (9) instead. We also consider this absorption in the 
following case studies. 

4. CASE STUDIES 

In this section we study the three bright Fermi-LAT GRBs 
0809 16C, 090510 and 090902B, and constrain their LFs with 
assumptions of one-zone or two-zone origins. 

4.1. GRB 08091 6C 

This is a bright long GRB, with a duration of ~ 50s and 
145 photons detected above 100 MeV, among which 15 are 
beyond 1 GeV and 1 beyond 10 GeV. The redshift is quite 
high, z = 4.35, so that the isotropic-equivalent energy is turned 
out to be E iso = 8.8 x 10 54 erg, the largest energy measured so 
far (Abdo et al. 2009a). 

The wide energy range spectrum of this GRB is well fit by 
a single Band function, which may imply that all radiation is 
originated from one region. Indeed with one-zone assump- 
tion, the one-component spectrum favors synchrotron origin 
over IC emission, and the spectral slopes can be understood 
in the frame work of synchrotron emission model (Wang et 
al. 2009). Using the time interval 3. 58-7. 68s, and under one- 
zone assumption, the LF has been constrained to be V > 900 
by Abdo et al. (2009a). It should be noted that the constraint 
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FIG. 3. — In the two-zone case the suppression factor /, due to attenuation 
by inner-originated and beamed target photons, as function of the observed 
photon energy e. Here the target photons are assumed to be a single power 
law distribution. All the lines correspond to the same parameters except those 
labelled in the plot. The dashed lines correspond to the totally beamed cases, 
calculated by using eq (10), while the solid lines are for the partly beamed 
cases using eq (12). 

is variability time dependent. In this constraint St = 2s is 
adopted from GBM light curve. However INTEGRAL also 
detected this GRB and show variability time in MeV range 
much shorter, St < 100 ms. With this shorter variability time, 
we constraint the LF to satisfy V > 1 130 (Fig 1). 

However, the onset of > 100 MeV emission is ~ 4 s de- 
layed relative to MeV emission; and in time bin "b" the bulk 
emission shifts toward later time as the photon energy in- 
creases, as pointed by (Abdo et al. 2009a), and the shift is 
1-s scale, much longer than the variability times in MeV light 
curves, St < 100 ms, as noted by Li (2010). These temporal 
behaviors suggest that HE emission may have different ori- 
gins and larger emission regions than the MeV one. Indeed, 
the "single" spectral component favors synchrotron over in- 
verse Compton radiation, however, as pointed by Li (2010), 
the observed highest energy photon in this GRB cannot be 
generated by synchrotron radiation, implying different com- 
ponent/origin for the HE emission. 

Here we consider a simple two-zone case, where the ejecta 
that produce HE emission > eq is released with a delay 
td = 1 - 4s relative to that produce MeV emission. It is hard 
to determine the threshold energy eq currently, but we take 
£o > 30 MeV, due to the different temporal behaviors above 
30 MeV. We use the observed flux and spectrum to calculate 
the optical depth due to absorption by inner-coming, beamed 
photons, and only use that at > eo to calculate the optical 
depth due to self absorption by local-originated photons from 
the HE emission region. With the sum of these two optical 
depths we can calculate the suppression / factor to modify the 
original HE emission that is free of absorption. We consider 
the time interval 3.58-7.68s following Abdo et al. (2009a), 
and assume the observed HE spectrum as the original one. 

The result is presented in Fig 4. It can be seen that, with 
r = 600 and eo = 30 MeV, the "self" absorption is less impor- 
tant than the "beamed" absorption, and the attenuated spec- 
trum does not show sharp cutoff but a slight steepening, as in 
Fig 3, in contrast with the one-zone case. We also show that 
taking eo ^ 30 MeV does not change the conclusion much. 
Indeed, for a photon of 3 GeV, the highest observed energy in 
the relevant time bin, and given T = 600 the threshold energy 



of 77 interaction is e th = T 2 (m e c 2 ) 2 /3GeV(l + z) 2 w 1 MeV, 
much smaller than 30 MeV. Moreover, we try different T val- 
ues and find the break energy, where the steepening happens, 
increases with T. T ~ 600 can be consistent with the observed 
spectrum in the two-zone case. Finally, it should be noted that 
the self-absorption becomes important when taking smaller 
threshold energy and time delay, i.e., eq = 10 MeV and t d = 1 s 
(for T = 600). 




e (keV) 

FIG. 4. — The HE suppression in the two-zone case of GRB 080916C. Here 
the observed, best-fit spectrum (see Table 1) is assumed to be the original 
spectrum without attenuation (dashed lines). The vertical line marks the ob- 
served highest photon energy. The calculated attenuations take into account 
the absorptions by both the inner- and local-originated photons. The upper 
panel shows the case of fixed threshold energy eo = 30 MeV, while the mid- 
dle and bottom panels show the attenuation varies with a range of threshold 
energy eo =10-50 MeV for delayed time of tj = 4 and 1 s respectively. In all 
these two zone cases V = 600 is taken. Note that the apparent dips in the last 
two panels show the contribution of self absorption. Also shown for compar- 
ison is the exponential cutoff (dot lines) in the one-zone case with T = 900 
and 8t = 2 s (Abdo et al. 2009a). 



4.2. GRB 090510 

This is a short GRB with a duration of 2. 1 s, but very bright, 
with 18 photons at > 1 GeV detected. Given the redshift z = 
0.903 ±0.003 and the total (0.5-1.0s) energy fluence in the 
10keV-30GeV band, (5.02 ±0.26) x 10" 5 erg cm" 2 , the total 
isotropic-equivalent energy release is (1.08 ± 0.06) x 10 53 erg 
(Abdo et al. 2009b). 

Using the spectrum in time interval 0.8 s-0.9 s which in- 
cludes a highest energy photon of 3 1 GeV and can be fitted 
by the Band function plus a power-law component, the LF 
constraint in one-zone case is T > 1200 (Abdo et al. 2009b). 
Under one-zone assumption we find that T > 990 (Fig 1). The 
two results are in broad consistence, though our result is a lit- 
tle less than that of Abdo et al. (2009b), which can be due 
to the different definition of R. Our defined R is larger by a 
factor of 2. 

Moreover, there are some distinct features in this GRB: the 
time-integrated spectrum in time interval 0.5- 1.0s is best fit 
by two spectral components, Band-function component at low 
energy plus power-law component dominating HE emission; 
the emission above 30MeV is delayed by td = 248 ms than 
those below 1 MeV as shown by the data analysis in Abdo 
et al. (2009b). These suggest that HE may have different ori- 
gin and/or emission region. Therefore we consider the simple 
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two-zone assumption for this GRB again. Since there are two 
components in the spectrum that may be consistent with the 
two components in the temporal behavior, we use the Band 
function component to calculate the optical depth due to MeV 
photon beam and use the power-law component for calcula- 
tion of self absorption. Thus we obtain the / factor to modify 
the power-law component, assuming that the observed best fit 
spectrum as the original one without 77 attenuation. 

The resulted spectra are shown in Fig 5. We can find that a 
LF of r ~ 600 can be still consistent with the observed spec- 
trum. It should be noted that although the power law com- 
ponent dominates in energy the Band function component 
still dominates in photon number. So the absorption due to 
beamed MeV photons can be more important. For the param- 
eters taken, the self absorption due to local originated photons 
contribute comparable, though less important, effect, thus the 
attenuated spectrum is steeper than beamed-Me V-photon only 
case, but still much smoother than the sharp cutoff in one-zone 
case. 
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FIG. 5. — The HE suppression in the two-zone case of GRB 090510. 
r = 600 is assumed, and t,t = 0.25 s is taken due to the data analysis by Abdo 
et al. (2009b). Here "self" denotes the absorption by local-originated photons 
in the HE emission region itself, "beamed" by inner-originated and hence 
beamed photons, and "beamed+self" by both inner- and local-originated pho- 
tons. Also shown for comparison is the one-zone case with T = 1200 and 
St = 12 ms (Abdo et al. 2009b). 



4.3. GRB090902B 

With the redshift of 1.822 this long, fairly strong GRB 
has an isotropic-equivalent energy Ei so = 3.63 ± 0.05 x 10 54 
erg, comparable with that of the highest-energy one GRB 
080916C (Abdo et al. 2009c). The duration in the energy 
interval 50-300 keV of Fermi (GBM) is 22 s. The highest 
energy photon (33.4 GeV) in this GRB is detected at 82 s af- 
ter trigger, while that in the prompt phase is 11.2 GeV and in 
interval of 9.6-13 s. Using this time interval and one-zone as- 
sumption, Abdo et al. (2009c) constrain the LF to be T > 1000 
(Abdo et al. 2009c), while we get, in Fig 1, V > 830. 

Similar to GRB 090510, this GRB also has a distinct spec- 
tral component fitted with a power law besides the Band func- 
tion one. A peculiar characteristic of its spectrum is that its 
power-law component extends to lower band (<10keV). Sim- 
ilar to GRB 0809 16C, there is an obvious delay of a few sec- 
onds in the HE onset. Look at the time bin "b" in Abdo et 
al. (2009c), the light curve peak seems also to shift toward 



high energy, with a one-second delay. Then we consider again 
a simple two-zone case taking td = 1 — 5 s. Similar to GRB 
090510, we use the Band function component as the beamed 
MeV target photons of two-zone absorption, and the power- 
law component for the self absorption at HE emission region. 
The total optical depth will lead to / factor calculation, which 
further modify the original spectrum, assumed to be the ob- 
served spectrum. The results in Fig 6 suggest that Y ~ 600 
can still be consistent with observations. 
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FIG. 6. — The HE suppression in the two-zone case of GRB 090902B. T = 
600 is assumed and the delayed time is taken to be tj = 5 s (upper panel) 
and 1 s (lower panel). Also shown for comparison is the one zone case with 
T = 1000 and St = 53 ms (Abdo et al. 2009c). 



5. DISCUSSION AND CONCLUSIONS 

We have revisited in this work the problem of constraining 
the GRB LFs by the HE attenuation. Although this problem 
has been considered by many previous works, two concerns 
that have been ignored in the previous work have been empha- 
sized here. First, we notice that in the one-zone case in order 
to self-consistently calculate the 77 optical depth one needs 
to consider the target photons with HE spectral cutoff, other 
than extending to infinity. This concern is important when 
the LFs are below a few hundreds, or when the luminosity of 
GRBs are low. Second, we relax the one-zone assumption and 
consider a simple two-zone case where the beaming of target 
photons in the emission region should be taken into account. 
Our results show that in the two-zone case, the 77 absorption 
does not lead to an abrupt spectral cutoff but a spectral steep- 
ening. If the target photon energy distribution is with a power 
law with photon index s then the spectral slope is changed by 
a factor of - — I. This also predicts that there should be no 
spectral cutoff in the GRB spectra if the prompt emission is 
not produced in one single region. 

It should be noted that there are some attempts by other 
authors to improve the approximation for the optical depth. 
Baring (2006) concluded that the pair attenuation signature 
appears as broken power-law rather than exponential cutoff 
by considering the skin effect and introducing an attenuation 
descriptor of 1/(1 +r) instead of e~ T . Granot et al. (2008) 
considered the emission zone as a very thin layer producing 
impulsive emission. They calculated in detail the opacity evo- 
lution during a pulse, and claimed that the attenuation signa- 
ture can be different from that derived from the simple one- 
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zone approximation. Essentially, these two works still con- 
cern one-zone problem, with AR ~ R Me y. However in the 
two-zone problem that we considered here, the HE and MeV 
emission components are emitted at very different radii, with 
AR 3> Rmzv, which leads to much smaller optical depth and 
hence smaller LF at HE emission region. 

Furthermore, we take our new concerns to analyze the spec- 
tra of the three bright GRBs 080916C, 090510 and 090902B 
and found that in the two-zone case a LF of V ~ 600 can still 
be consistent with the observed spectra. This relaxes the strict 
requirement, T > 10 3 , in one-zone assumption. 

We note that in the present observational situation where 
only tens to hundreds HE photons detected in one GRB, a 
slight change of the spectral slope is not easy to be identified. 
A single power law may still fit the HE spectral tail. 

We have considered a simple two-zone case here. However 
the situation can be more complicated. The central engines 
of GRBs may naturally create variabilities in a wide range of 
timescales, e.g., from ~ 1 ms to <~ 10 s. In the framework 
of internal shock model, this will lead to kinetic dissipation 
in a wide range of radii. Even in the single-timescale case, 
the internal collisions will happen as the ejecta expand until 
the material is distributed with velocity increasing with ra- 
dius. In such case we will expect multi-zone other than sim- 
ple two-zone case. The time-integrated spectrum- note that 
the time interval with high enough photon statistic is usually 
much larger than the variability time- will be contributed by 
the multiple regions. We also calculate cases with T < 600, 
the sum of the flux at the HE end can be comparable to the 
original flux. This means that the spectra can be consistent 
with a multi-zone case with the LF T < 600. 

The formation and acceleration of relativistic collimated 
GRB jets are open questions. In the standard "fireball" model, 
the thermal pressure can only accelerate the gas up to a LF 
T < 10 3 (see, e.g., Piran 1999; Li 2010). On the other hand, 
simulations of magnetic-driven jets (e.g., Tchekhovskoy et 
al. 2009) can generate jets with the product of the LF and 
jet opening angle being T0j « 10-30, which is consistent 
with pre-Fermi GRB observations. However for the bright 



Fermi-LAT GRBs, Cenko et al. (2010) constrained the jet 
opening angles by their afterglows, which, combined with the 
large LF, T > 1000 from 77 attenuation argument, suggests 
much larger values of TOj. We stress here that if relaxing the 
one-zone assumption for GRB multi-band emission, LFs with 
"normal" values, say, T < 600, can still be consistent with ob- 
servations. This relaxes further the theoretic problem of jet 
acceleration. 

Recently, a similar paper, Zou et al. (2010), considering 
the same two-zone absorption, is now in preprint. The main 
difference between two papers is the rest frames for the op- 
tical depth calculation, i.e., we consider the comoving frame 
of the HE emission region while they consider the observer 
frame. They integrate over the region up to R mdx , where the 
HE photon spatially leaves the MeV front, while our integra- 
tion corresponding to one dynamical time expansion is equiv- 
alent to integration up to 2R, where the generated HE photon 
doubles its radius. Because both the number density of the 
target photons and the angle between the travelling directions 
of HE photon and the MeV front decrease rapidly with radius, 
the interaction is strongly dominated by those at small radius, 
and hence the upper limit of the integration is unimportant- no 
matter the upper limit is R mix ('^- R) or 2R the result is prac- 
tically the same. Furthermore, they use an "averaged" optical 
depth to constrain the LF, which may not be appropriate since 
we have shown that no sharp cutoff is expected in the two- 
zone case. We consider more carefully the spectral profile 
due to suppression. Finally, they neglect the self absorption in 
the local region which may contribute significant effect as we 
show. 
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APPENDIX 
DERIVATION OF THE SUPPRESSION FACTOR FOR THE TOTALLY BEAMED CASE 

Here we derive the suppression factor / (eq. 19) in the two-zone case, assuming the target photon distribution as a single power 
law with N(e) = Noe~ s . In the comoving frame of the HE emission region, the target photon distribution can be given by eq (4). 

The cross section of 77 collisions in the relativistic limit (E ^> m e c 2 ) is a(E) w <7qE~ 2 , where cto = 
(3/8)crr(m e c 2 ) 2 [21n(2£ , /m e c 2 )- 1] weakly depend on E and can be considered as constant due to roughly a constant of 
2ln(2E /m e c 2 )- 1, which is as an approximation taken as 3. With these approximations the 77 optical depth (eq.12) can be 
reduced to 



t = CiT*- 3 e"-\l - m')' ? (1 +zf-\ C x = Noa Q (2m 2 e c 4 r s d 2 L (sc 2 t d T l . 
Using the transformations of eqs (13) and (14), the optical depth further becomes 

r = C 2 eT\l - M ) s (l -/SrM)"' , C 2 = Ci2'(l + zf^T 4 . 
As 9 increases (/1 decreases) r increases. Let us define the minimum /i m ; n where r(e,^ m i n ) = 1, then at \i < [i 



(Al) 



(A2) 



at e is significantly absorbed. Using the approximation that e T 
factor is approximated by 



min the emission 
1 when t < 1 and e~ T = when r > 1, the expression for / 



f(e;T)^f dutl-friM) 



-h-\ 



(l-^rr-'a-Mmin). 



The second equality holds for /i m ; n — ► 1 . 
From the definition of /i m i n , we can solve out /x n 



1 -Mmin = C 2 ' (1 - /3rMmin) 7 £ 7 ~ Q> ' (1 - /3 r ) 7 



1 + 



/3r(l-Mmin) 



s(l-/3r) 



(A3) 



(A4) 
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where the second equality in the above equation, again, holds for /i m ; n — *■ 1 . 
We can see that approximately, 

f - 1 " 1 l-Mmin<^(l-/3r) 



A /^min ^ 



const. 1-Mmin > -sril-fir) 



■k 



The critical condition 



l-Mmin=-5-(l-/3r) 

PT 



(A5) 



(A6) 



corresponds to a "break energy" e\, v . Substituting eq (A6) into the both sides of the first equality in eq (A4), we get 



£br 



s(l +zY 



(i+s)KcT4 



N a dl 



e > £br corresponds to 1 -^ m ; n < -£-(l—fir) and e < £br to 1 -^ m ; n > -S-(l-(3r)- So we can write 



f(s;T)(X l-/i min oc 



[const, e < £b r 



or, after an arbitrary normalization, 



/(£;i> 



£ > £br 



\ e *<) 



£< £br 
£ > £br 



(A7) 



(A8) 
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